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Abstract. We consider a modulated discrete nonlinear Schrodinger (DNLS) model with alternating on-site 
potential, having a linear spectrum with two branches separated by a 'forbidden' gap. Nonlinear localized 
time-periodic solutions with frequencies in the gap and near the gap - discrete gap and out-gap breathers 
(DGBs and DOGBs) - are investigated. Their linear stability is studied varying the system parameters 
from the continuous to the anti-continuous limit, and different types of oscillatory and real instabilities 
are revealed. It is shown, that generally DGBs in infinite modulated DNLS chains with hard (soft) non- 
linearity do not possess any oscillatory instabilities for breather frequencies in the lower (upper) half of 
the gap. Regimes of 'exchange of stability' between symmetric and antisymmetric DGBs are observed, 
where an increased breather mobility is expected. The transformation from DGBs to DOGBs when the 
breather frequency enters the linear spectrum is studied, and the general bifurcation picture for DOGBs 
with tails of different wave numbers is described. Close to the anti-continuous limit, the localized linear 
eigenmodes and their corresponding eigenfrequencies are calculated analytically for several gap/out-gap 
breather configurations, yielding explicit proof of their linear stability or instability close to this limit. 

PACS. 63.20.Pw Localized modes - 63.20.Ry Anharmonic lattice modes - 42.65 Nonlinear waveguides 



1 Introduction 

The existence of spatially localized nonlinear excitations 
in discrete models - 'discrete breathers' (DBs) - has at- 
tracted much attention, as they can play a significant role 
in condensed matter physics, biophysics, nonlinear optics, 
etc. (for reviews see e.g. mH]). In particular, being local- 
ized in space and typically stable, DBs can contribute to 
energy and information transfer processes. For example, 
breathers can trap energy during long periods of time and 
cause non-exponential thermal relaxation in nonlinear lat- 
tices 0. By contrast, under certain conditions 4t|5 DBs 
can be rather mobile and become good energy carriers. 

The presence of intrinsic structure of a medium can in- 
fluence the DB properties. In particular, new types of DBs 
- discrete gap breathers (DGBs) - appear with frequencies 
inside the forbidden gaps in the linear waves spectrum. 
DGBs are the discrete analogues of gap solitons, first dis- 
covered in a nonlinear optical medium with modulated 
refractive index 6 . Their specific features are due to the 
existence of two neighboring bands of the linear dispersion 
curve, with opposite signs of dispersion close to a gap. 

A classical example of a system with two bands in 
the linear spectrum is a diatomic one-dimensional lattice. 
DGBs in diatomic lattices with nonlinear interatomic po- 
tentials [Fermi-Pasta-Ulam (FPU) models] were studied 
analytically as well as numerically within the rotating- 



wave approximation, neglecting higher harmonics gener- 
ation [7|IH|F?)ll()llll| . Later the existence of DGBs in di- 
atomic FPU lattices was proved rigorously 12.13) . Exis- 
tence and linear stability properties of DGBs in diatomic 
FPU chains were also studied numerically for some partic- 
ular values of system parameters (14ttl5al6) . We recently 
[T7j performed a detailed analysis of DGBs in a diatomic 
lattice with nonlinear on-site potential - Klein- Gordon 
(KG) model - in the complete regime of continuation from 
the anti-continuous (AC) to the continuous limit, and in 
particular we described the dynamics resulting from sev- 
eral types of instability mechanisms (oscillatory as well as 
non-oscillatory) . 

Very recently, much attention has been attracted to 
modulated structures described by equations of discrete 
nonlinear Schrodinger (DNLS) type, in particular within 
nonlinear optics for the description of electromagnetic waves 
in arrays of weakly coupled optical waveguides '18'. For 
such a system with even and odd waveguides of different 
widths, discrete gap [T^ll^ as well as multigap soli- 
tons were found. A related model with two coupled DNLS- 
like equations describing a diffraction-managed waveguide 
array was also studied recently, and found to support dif- 
ferent types of discrete gap solitons j22j . Much recent at- 
tention has also been given to DNLS models in describing 
Bose-Einstein condensates (BECs) trapped in optical pe- 
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riodic potentials |23II24| . For a recent review on properties 
and applications of the DNLS equations, see |25j . 

As the DNLS equations also approximate the small- 
amplitude dynamics of weakly coupled KG-chains [261127) , 
the DNLS discrete gap solitons should share many features 
of DGBs in KG (and also FPU) lattices (henceforth they 
will therefore also be referred to as DGBs). In particular, 
in the continuum limit these models are described by the 
same general class of coupled equations with exact gap 
soliton solutions analyzed e.g. in |28II29| . 

The aim of the present study is to give a detailed and 
thorough analysis of DGB properties in the modulated 
DNLS model. Our motivation is twofold. First, we wish 
to provide, for the benefit e.g. of experimentalists work- 
ing in related areas (nonlinear optics, BECs), charts of 
different types of instabilities etc., and information about 
the parameter regimes where they are expected to be seen. 
Particularly interesting is the possible observation of mov- 
ing DGBs. Generally DBs can not freely propagate along 
the system due to the existence of the effective periodic 
Peierls-Nabarro (PN) potential (caused by the discrete- 
ness of the system). However, in the diatomic KG model 
we found ^Tj parameter regimes with increased DGB mo- 
bility due to 'exchange of stability' between symmetric 
and antisymmetric DGB configurations, and thus it is 
of interest to look for similar effects in the modulated 
DNLS models, proposing that the increased DGB mobility 
should be directly observable in experiments. 

Second, it was not carefully analyzed before, within a 
fully discrete model, what happens with DGBs when their 
frequencies approach gap boundaries. In continuous mod- 
els, gap solitons delocalize and vanish at one gap bound- 
ary, but bifurcate into a new type of excitations - out- 
gap solitons - at the other boundary ^Hl- Out-gap soli- 
tons have frequencies inside the linear spectrum, and their 
structure can be viewed as a superposition of two fields. 
One field has a 'dark' soliton shape with non-zero ampli- 
tudes at the infinities and a localized decrease in the cen- 
ter. The other field has either a 'bright' soliton form with 
tails exponentially decaying to zero, or forms a 'bright soli- 
ton on a pedestal' ('anti-dark soliton'), similar to a 'dark' 
soliton but with a localized increase of the amplitude in 
the center. Although out-gap solitons can radiate energy 
through linear waves, they may still represent long-lived 
localized excitations persistent to perturbations (e.g. j31p. 

Thus, analogously it is interesting to look for discrete 
out-gap breathers (DOGBs) in discrete models, and to 
study their stability and the bifurcations from DGBs to 
DOGBs. This was not analyzed in earlier studies of DGBs 
in KG or FPU models, partly due to technical problems: 
since DGBs delocalize when approaching gap boundaries, 
numerical computations become time consuming when in- 
creasing the system size to avoid boundary effects. How- 
ever, DNLS models are more conveniently analyzed nu- 
merically, since time-periodic DGBs become 'stationary' 
solutions with purely harmonic time-dependence, that can 
be removed by transforming into a rotating frame. 

Another important question concerns the possible exis- 
tence of truly localized excitations with frequencies inside 



the linear wave spectrum - analogues of 'embedded' soli- 
tons (e.g. JS2j) - and the nature of tails in such hypotheti- 
cal localized DOGBs. This issue is also most conveniently 
analyzed with the modulated DNLS model. 

Our study of the spatially binary modulated DNLS 
model and the properties of its DGB and DOGB solutions 
will be structured as follows. In Sec. [3 the model is de- 
scribed, as well as the properties of its linear spectrum. In 
Sec. 13 we discuss the linear stability analysis of stationary 
solutions, and in Sec.0]we review the numerical construc- 
tion of breather solutions from the anti-continuous limit, 
discussing some particular features of the continuation in 
couphng and frequency of DGBs and DOGBs. In SecElwe 
present numerical and analytical results for linear stability 
properties of different types of DGBs and DOGBs. Sec. 15. II 
describes the continuation in coupling of DGBs, and we 
compare the stability results with those obtained for the 
diatomic KG chain TT. In Sec. 15.21 wc study the transi- 
tion from DGBs into DOGBs while changing the breather 
frequency, describing the bifurcations which may occur to 
DGBs when their frequency approach the gap boundary, 
and furthermore penetrate the linear spectrum as DOGBs. 
In Sec. 15. HI bifurcations of DOGBs with different types of 
tails are described. We study the subsequent bifurcations 
of DOGBs into 'on-top' DBs (DOTBs) with frequencies 
above the linear spectrum, and obtain a general bifurca- 
tion picture of DOGBs with tails of different wave num- 
bers. We discuss the linear stability properties of the most 
important DOGB and DOTB solutions in Sec. 15.41 and in 
Sec. El conclusions are made. Some details of the analyti- 
cal investigation of linear stability properties of DGBs and 
DOGBs at small coupling are deferred to Appendix El 



2 The Model 

The standard DNLS equation with cubic nonlinearity reads: 

= ^V'n - C {lpn-l + V'n+l " 2?/^„) -|- 7 \^nf (1) 

When modelling an array of weakly coupled optical wave- 
guides |18II19| , V'n is the normalized amplitude of the elec- 
tric field in the EM wave, propagating along the array. The 
coefficient A then characterizes the linear propagation con- 
stant in the waveguide, C the coupling coefficient, 7 the 
effective nonlinear coefficient, and the time-like variable t 
measures distance along the array. Here we consider pos- 
itive 7, corresponding to 'hard' nonlinearity. The case of 
'soft' nonlinearity (7 < 0) can easily be recovered by a 
simple transformation, see footnote |H below. 

To introduce a gap parameter, one can modulate any 
of the constants A, C, 7 (or two or three of them simulta- 
neously). We choose to modulate the coefficient A, corre- 
sponding in the optical model to a variation of the width 
of even and odd waveguides (see e.g. [I9|). Thus, we put: 

X-^ Xn = l + {-irS\ (2) 

where the parameter S defines the width A of the gap in 
the linear waves spectrum: A = 26"^ (see below). 
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Fig. 1. The dispersion relation of linear waves 

The dispersion relation for solutions: 

. I Aexp{iqn ^ iutt), n^2k, , . 

\Bexp{iqn-iujt), n = 2k + l, 

fc = 0,l,2,... 

of the linearized system ^ with A in the form Q reads: 
LUo,a = 1 + 2C± ^S^ + 4C^ cos2(q), (4) 

where indices o and a stand for the upper ('optic-like') 
and lower ('acoustic- like') branches of the spectrum^. It 
possesses a gap at the wave number q = n/2 (see Fig.P). 
The gap boundaries are defined as: 



CJ2 1 = 1 + 2C ± 



(5) 



The ratio of the amplitudes A and i? in a linear wave ^ 
is defined by its frequency uj: 



A 
B 



(6) 



In linear waves with frequencies belonging to the upper 
branch of the spectrum (lo > UJ2), amplitudes in even 
waveguides {A) are always higher than in odd waveguides 
(B), while for frequencies in the lower branch (w < wi) 
the odd waveguides are dominant. Therefore, the sub-field 
of even (odd) waveguide amplitudes can be referred to as 
the upper (lower) band sub-field, respectively. With an 
analogy to a diatomic chain, even and odd waveguides 
correspond to light and heavy atoms, and will therefore 
be referred to as 'light' and 'heavy' sites, respectively. 

Another type of modulated DNLS model is obtained 
by multiplying the left-hand side of equation with to„ , 
defined analogously to A„ in (0) . This model, correspond- 
ing to simultaneous modulation of all constants in , can 
be considered as a simplified model of a diatomic lattice 

^ Note however that also for the "acoustic-like" branch the 
dispersion is quadratic and not linear for small q. 



(rrin playing the role of atom masses), neglecting higher 
harmonics generation. As shown in for the diatomic 
FPU model, the effect of higher harmonics on DB prop- 
erties is often rather negligible, although additional insta- 
bilities and bifurcations generally may appear for large 
values of the coupling and/or large- amplitude oscillations 
(e.g. 27'). Investigating numerically breather properties in 
the two modulated DNLS models we have found qualita- 
tively similar results, and therefore we focus our discussion 
only on model with coefficients A as defined in 

We denote the system size (the number of waveguides) 
with N {n = 1,2, N) and impose periodic boundary 
conditions ipN+i = V'lj "00 = V'w- The DNLS equations 
are then the Hamiltonian equations with the Hamiltonian: 



^^({#4,{C}) 



(7) 



N 



J2 (C|V„ + l-V'n|'+A|V„|' + ||V'n| 



In what follows we put 7 = -fl in (Q, (TJ (without loss of 
generality, since varying 7 is equivalent to rescaling ipn)- 

3 Stationary solutions and linear stability 

Stationary solutions of equations 1^ have the form: 

^Pn{t) ^ (j)ne-'^'\ (8) 

with time-independent amplitudes (j)n , where cob is the ex- 
citation frequency. DGBs have frequencies Ub lying inside 
the gap (O, so that ui < ojb < uj2- For DOGBs one has 
u!b > 0J2, while 'on-top' breathers (DOTBs) have frequen- 
cies above the linear spectrum, LUb > where 



w„ = 1 + 2C + VS^+iC^. (9) 

To analyze the linear stability of a particular solution 
{</'«'''}, we add a small perturbation {e„(t)} to it: 

Mt) = {cb^^'>+en{t))e-^-^\ (10) 

and linearize the equations 

= + - ^"-C (e„_i + e„+i) + (11) 

at 



+27 



/,(o) 



Expanding e„(i) in real and imaginary parts, e„(t) = a„(t)- 

iPn{t), equations l(TT|l can, for real solutions {(jin^}, be 
written in matrix form |34| : 

'{a„}\ _( Co\( 

{/?„>; y-c, J 

with £0 a-nd Ci defined as follows: 



(12) 



\.a + 2C-UJb+l(dp^n^)^ 

Xn + 2C-LJb + 3j (4"^) 
-C (a„-i + a„+i) . 



f3n 



(13) 
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The solution {(f>n } is linearly stable if and only if the per- 
turbation {(-nit)} remains bounded in time, i.e. if and only 
if all eigenvalues iuj^ of the matrix in (|12|l are imaginary. 
Note that writing the general solution of as 

= i («n + bn) e--* + li<- K) e''^'' (14) 
leads to the following set of algebraic equations (e.g. 27j: 

\ £ia„ = OJebn, ^ ' 

with £o and Ci defined above l(T^. Consequently, the con- 
dition of the linear stability of a solution {(^n^"*} is that all 
the eigenfrequencies Wg, obtained from H15|) . are real. 

If We is a real eigenfrequency of (|15|l and (a„, 6„) is the 
corresponding real solution, then ioje is an eigenvalue of 
ifT^ with eigenvector ({a„ = a„}, {/?„ = i5„}). Therefore 
with each real eigenfrequency We > one can associate a 
Krein signature ^(we) (e.g. 35:,27^ and references therein): 

nioje) = sign^ a„&„. (16) 

n 

For DGBs with frequencies inside the gap, eigenmodes 
with positive (negative) n correspond to excitations inside 
even (odd) waveguides in the AC limit C ^ 0. Physically, 
the Krein signature is the sign of the Hamiltonian energy 
carried by the corresponding eigenmode (e.g. jSSlEni)- 

4 Construction of breather solutions 

The main idea for the numerical methods used to obtain 
breather solutions, based on the proof of their existence 
[Tll37j ■ is to perform numerical continuation of an exact 
breather solution, known for some particular parameter 
values, to other, arbitrary, values of these parameters. For 
stationary solutions (jHJ, Eq. |^ with modulated coeffi- 
cient A (0) reduces to a set of algebraic equations (7 = 1)^: 

[l + (-l)"(52-cc>,]0„- (17) 

-C (0„-l + (f)n+l - 20„) -t- \(j)n\^ 4>n = 0, 

and solutions can be found by Newton schemes (e.g. 38 ). 

Usually breather solutions at C = (AC limit) are 
taken as initial solutions and continued to non-zero values 
of C. At C = one has an array of uncoupled waveguides. 
Hence, in the AC limit, for real solutions each amplitude 
0„ should be taken from the set (/)„ £ {0,±yl„}, with: 

An = y/uJb-Xn- (18) 

Consequently, to each breather solution one can associate 
a corresponding coding sequence, where we use symbol 

^ Note that the case of soft nonlinearity 7 = — 1 is also 
covered through the transformations (j}'„ = ( — oj^ — 
2 4- 4C — ujt,, 5'^ = —5^. Essentially, this reverses the roles 
of upper/lower bands and even/odd sites. 



notation with up- and down- arrows for +An and —An, 
respectively. To make the intrinsic structure of the solu- 
tions more clear, we use different sets of symbols for odd 
and even amplitudes (lower and upper band sub- fields): 
{O trJI} and {0 ti}, respectively. 

Note, however, that these coding sequences are not 
unique for a given solution. Indeed, one can construct a 
DGB with some particular coding sequence and then con- 
tinue it in frequency up to values above the upper gap 
boundary Ub > UJ2- Above the gap this solution will have 
another coding sequence if continued to C = 0, though 
the continuations in frequency and coupling can be done 
smoothly without any true bifurcations to other solutions. 
Similar effects were observed for second-harmonic pene- 
tration into a phonon band for breathers in a monoatomic 
Klein-Gordon lattice [HHI- In fact, one should distinguish 
two types of coding sequences: those corresponding to the 
AC limit inside and above the gap of the linear spectrum 
('gap' and 'out-gap' coding sequences, respectively)'^. We 
will use the superscripts G or O in coding sequences to 
indicate whether they are 'gap' or 'out-gap', respectively. 

Restricting to solutions with spatial symmetry or anti- 
symmetry with respect to the center, Nc = Int(A^/2) -I- 1, 
the coding sequences can be simplified by omitting the 
codes for sites n — l,2..Nc — 1, and using subscripts S 
or A to mark whether a breather has symmetric or an- 
tisymmetric spatial configuration respectively. As follows 
from 1)18(1 and (0) , in the case of hard nonlinearity (7 > 
in (Q), the simplest symmetric DGBs are centered on a 
'heavy' (odd) site with non-zero amplitude, while the sim- 
plest antisymmetric DGBs are centered on a 'light' (even) 
site with zero amplitude and with non-zero amplitudes of 
opposite signs on two neighboring 'heavy' sites. 

Note that the structure of the linear spectrum, and 
in particular the gap boundaries uji^2 ®; depends on the 
coupling C (although the gap width A = 26"^ is indepen- 
dent on C). Thus, if one constructs a DGB with frequency 
ojb inside the gap at some particular value of C and tries 
to continue it versus C at fixed frequency, it can reach 
one of the gap boundaries in the process of continuation. 
Consequently, a DGB can disappear (on the lower bound- 
ary) or transform into a DOGB (on the upper boundary) 
while continued in coupling. The same effect can occur 
for DOGBs, as the frequency may reach that of the linear 
wave corresponding to its tail. To avoid this problem we 
fix not the breather frequency, but the frequency detuning: 

AtJ = uJb-uJo{q,C), (19) 

for the continuation in coupling. Here tOoiq, C) is the fre- 
quency of the linear wave Q with the same wave number q 
as that of the breather tail, belonging to the upper branch 
of the spectrum. All stationary DGB solutions have tails 
with wave number q = 7r/2, so that the frequency detuning 
(tl^ is equal to Aoj = oJb — 0J2{C) < for such breathers. 

^ The existence of two types of coding sequences is connected 
to the existence of two types of AC limits. In general, the AC 
limit uJb/C 00 can be reached either by taking C = or 
LJ6 — > 00. These are not equivalent for modulated systems. 
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Hence, the frequency detuning Auj (|19(l is a more con- 
venient parameter than the frequency LUf,, and it can be 
considered as the only dynamical parameter of station- 
ary DGBs and DOGBs. Note also that although Eq. IH} 
for stationary solutions has three parameters (wb, S^, C), 
only two are independent, since varying the third param- 
eter is equivalent to rescaling ipn- For instance, one can 
choose as independent parameters: S'^/C - the 'modula- 
tion' parameter, and uJb/C - the 'discreteness' parameter. 
In what follows we fix the 'gap' parameter 5^ and inves- 
tigate the properties of stationary solutions while varying 
the frequency detuning Au) H19|l and the coupling C. 



5 Numerical results 

We here present the results of numerical continuation and 
stability analysis of DOB and DOGB solutions. 

5.1 Stability of gap breathers 

At frequencies inside the gap, oji < cob < ^2, only 'heavy' 
(odd) amplitudes (/)„ can be non-zero in the AC limit (|18|) . 
Thus, there are two basic configurations of DGB solutions 
with coding sequences {i[ (0O)}f and {0 i[ (OO)}^ (codes 
in parenthesis are repeated), representing a symmetric 
DGB centered on a 'heavy' site with non-zero amplitude, 
and an antisymmetric DGB centered on a 'light' site with 
zero amplitude, respectively (Fig.|3). Since DGBs have ex- 
ponentially decaying tails of wave number q = n/2, neigh- 
boring even amplitudes, as well as neighboring odd am- 
plitudes, have opposite signs. In addition, there is a phase 
shift of 'light' (even) amplitudes in the center (Fig. 
known from gap solitons in the continuum limit (e.g. 29 ). 

In Fig. 121 we show the results of stability analysis of 
symmetric and antisymmetric DGBs for the continuation 
in coupling at three different frequency detunings (fT^ : 
in the upper half of the gap, in the gap center, and in 
the lower half of the gap, respectively. The stability eigen- 
values can be divided into two groups, corresponding to 
spatially localized and extended eigenvectors, respectively. 

The 'extended' eigenvalues are independent on the brea- 
ther spatial configuration and would be obtained also for 
the system without breather. They correspond to extended 
excitations of the upper and lower band sub-fields and 
form two bands icu^ and iw- with positive and negative 
Krein signatures, respectively: 

iuj^ = i (ujoiq) - ujb) , (20) 
ztj_ ^ i{ujb~ ujaiq)) ■ 

Here Wo((?) and uja{q) are the frequencies of linear waves 
lO, O with wavenumbers q (defined by the system size), 
corresponding to the upper and lower bands, respectively. 
In the limit of an infinite system the extended eigenvalues 
will have continuous spectrum. The boundaries of these 
bands are determined by the DGB frequency ujb and the 
boundaries of the lower and upper bands of the linear 
spectrum u>i^2,^u.d (see Fig. Defining the deviation 



a = Auj + of a, DGB frequency from the gap center: 
tJ6 = 1 + 2C + a {-S^ <a< S^), yields 

~a<uj+<VS'^ + 4C2 - a, (21) 

+ a. 

Thus, having the same shape, the two bands of ex- 
tended eigenvalues are shifted by 2a with respect to each 
other. While increasing the coupling constant C, these 
bands broaden and at a certain value Co, 

Co = V^^RT^, (22) 

they overlap. (In the gap center a — 0, the two bands co- 
incide as shown in Fig. intc),(d).) The overlapping of the 
bands will lead to collisions of extended eigenvalues with 
opposite Krein signatures with generation of instabilities 
(see Fig.OJ. For larger systems more collisions will occur, 
as the bands become more 'dense', but the instabilities 
produced will become weaker and should completely dis- 
appear in the limit of an infinite system 40 

The localized eigenvalues play a more important role, 
as their collisions with each other or with extended eigen- 
values can produce strong instabilities independent on the 
system size. The corresponding eigenvectors are localized 
around the breather when the eigenvalues lie outside the 
extended bands. There are two types of localized eigenval- 
ues, bifurcating from the top of the extended band with 
K = +l (PI, P2, P3, ... in Fig.EIa)), and from the bottom 
of the band with k = —1 {Nl, N2), respectively. In addi- 
tion, there is also the 'phase' mode with itUe = 0, whose 
eigenvector is the same as the breather solution, and as- 
sociated with the phase rotation of the breather. 

The eigenvector corresponding to the eigenvalue A^l 
has symmetry opposite to the DGB. Colliding with its 
complex conjugate a real instability is produced, connected 
to the 'exchange of stability' between symmetric and an- 
tisymmetric DGBs ^ [Peierls-Nabarro (FN) mode, see be- 
low]. The other eigenvalue bifurcating from the negative 
Krein signature band, N2, has an eigenvector with same 
symmetry as the breather. For DGB frequencies in the up- 
per half of the gap, both iVl and N2 eigenvalues penetrate 
the positive Krein signature band, yielding oscillatory in- 
stability (Fig.|3{a), and 7V2 in (b)). In particular regimes 
the N2 eigenmode can collide with its complex conjugate 
producing real instabilities, connected to transitions from 
'discrete-like' to 'continuous-like' DGBs jl7j (see Sec. 15. 21 ) 

The eigenvalues P1,P2,P3, ... bifurcate one by one 
when increasing the coupling, corresponding alternately 
to symmetric or antisymmetric eigenvectors. The number 
of these localized modes increases as, with increasing cou- 
pling, the central part of the DGB will occupy a larger 
number of sites. These eigenvalues can also produce os- 
cillatory instabilities in the upper half of the gap, pene- 
trating the extended band with opposite Krein signature 

* This is true generically for localized breathers, while for 
non-localized solutions (e.g. breathers with extended tails) 
such instabilities survive also in the infinite-size limit (e.g. 127]). 

^ Here 'stability' only means w.r.t. the PN-mode, since Krein 
instabilities may exist also in this regime (see Fig. |3a),(b)). 
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Fig. 2. Symmetric (a) and antisymmetric (b) DGBs at non-zero 

(see Fig. |3Ia),(b)). Unlike other eigenvalues, PI and P2 
(and for the antisymmetric DGB also P3) can bifurcate 
from the extended band immediately when C becomes 
non-zero. Analysis shows, that for each of these eigenval- 
ues there is a critical value of the breather frequency de- 
tuning AiOcr (in the lower half of the gap for PI and P2), 
specific for symmetric and antisymmetric DGBs, above 
which they stay outside the band (i.e., the eigenmodes 
are localized) at arbitrarily small non-zero coupling (see 
Appendix 1X1 for details). The deviation of these eigenval- 
ues from the top of the extended band is proportional to 
at small values of C (see Eas.(IS ^ . (|S5|) '). By contrast, 
below Aixicr the eigenvalue is inside the band at small C . 

Similar structure of extended and localized eigenvalues 
was reported for the diatomic KG model ^17.. However, 
an important difference is that in the modulated DNLS 
model, DGBs with frequencies in the lower half of the gap 
do not possess any oscillatory instability associated with 
localized eigenvalues. Indeed, for frequencies in the lower 
half of the gap (a < in (|21|) '). the bottom of the extended 
band with n = —1 is below the bottom of the k = -1-1 band 
for all values of C (see Fig.|3Ic)-(f)), so the A''2 eigenvalue 
does not penetrate the extended band. For similar reason 
the localized eigenvalues P1,P2,P3... do not penetrate 
the extended band, as the top of the n = +1 band is above 
the top of the other band for all C. By contrast, in the 
diatomic KG model DGBs possess oscillatory instabilities 
associated with localized eigenmodes of P-type both in 
the upper and lower halves of the gap for larger values of 
the coupling, as the bands have different shapes |17j . 

The two types of oscillatory instabilities of DGB solu- 
tions were recently observed also in [20] , where the stabil- 
ity of DGBs was analyzed for a continuation in frequency 
at a particular value of coupling within a similar model 
of a coupled waveguide array. The P1,P2,P3,... insta- 
bilities can be associated with external resonances^ as the 
frequencies of the corresponding localized modes are above 
the linear spectrum. By contrast, the A'^2 instability can 
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coupling. Circles (triangles) denote 'heavy' ('light') amplitudes. 

be associated with internal resonances, as the frequency 
of this mode is inside the gap PU] , 

The real instability A^l was also observed in [201 for 
the antisymmetric DGB solution, which was shown to be 
unstable at all frequencies inside the gap for the particular 
used value of coupling. However, we emphasize that in fact 
both symmetric and antisymmetric DGBs can possess this 
type of instability (see Fig. at different values of cou- 
pling. This exchange of stability between symmetric and 
antisymmetric modes is very important for breather mo- 
bility issues |^17j . Generally, for the antisymmetric DGB 
at small C > there are two localized modes associated 
with antisymmetric and symmetric coupled oscillations in 
the two 'heavy' sites with non-zero codes. The former is 
the 'phase' mode with iuje = 0, while the latter is the A^l 
mode which becomes unstable immediately when C > 
(see Fig. |SIb),(d),(f) and Appendix IXjl . Thus, antisym- 
metric DGBs always possess a real instability at small C, 
with a growth rate proportional to C that can be obtained 
analytically f Appendix IXlTablelTl eigenmode (iVl, A, G)). 



5.2 Transition from DGB to DOGB solutions. 
Discrete- and continuous-like breathers. 

To study in details the transformation of DGBs into DOGBs, 
we perform the continuation in the frequency detuning 
(|19|l of the breather solutions, at fixed coupling. We focus 
our attention on symmetric breather configurations. 

In Fig. 01 the profile of a symmetric breather is shown 
for the continuation in frequency detuning, at fixed cou- 
pling, from a value inside the gap to a value above the lin- 
ear spectrum. At the upper gap boundary {Auj = 0) the 
DGB transforms into a DOGB with non-zero tails. The 
structure of this DOGB is qualitatively the same as that of 
a DGB (see Fig.l^a)) but with non-zero asymptotes of the 
upper band sub-field ('hght' amplitudes) |^^"^| ~ VAuj, 
while the lower band sub- field ('heavy' amplitudes) still 
exponentially decays to zero. The wave number of the tails 
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Fig. 3. Imaginary parts of eigenvalues iive of 11211 for symmetric (a),(c),(e) and antisymmetric (b),(d),(f) breathers at different 
values of the frequency detuning 11911 : (a),(b) Auj — —0.75 • (upper half of the gap); (c),(d) Auj — —S^ (middle of the gap); 
(e),(f) Au = -1.25 • 5^ (lower half of the gap). Light- grey crosses correspond to eigenvalues with positive Krein signature, 
dark-grey triangles to eigenvalues with negative Krein signature, and black stars to unstable eigenvalues. In all cases the system 
size N = 242 and the gap parameter 5^ = 0.2. 
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Fig. 4. Amplitudes of the central part of a breather for the con- 
tinuation in frequency detuning Au) from DGB [Auj = 
to DOGB {Auj = 5^). N = 50, 5^ = 0.1 and C = 0.02. 



H 




Aco 



Fig. 5. The dependencies of the Hamiltonian energy Q of 
a symmetric breather on its frequency detuning Au 1191 1 for 
different values of the coupling (from top to bottom): C — 0.17, 
C = 0.15 and C = 0.1. In aU cases = 0.1 and iV = 242. 

in this DOGB is q = 7r/2 as for a DGB, and the corre- 
sponding coding sequence is {fT (T O i In what 
follows we will call this configuration a 'normal' DOGB. 

In Fig. 13 we show the dependencies of Hamiltonian 
H ^ on frequency detuning Auj I|19() . for continuation 
of symmetric breathers from Auj = —26^ (bottom of the 
gap) to Auj = 0.1(5^ (inside the upper band) at different 
fixed values of the coupling. At the upper gap boundary 
Auj = 0, a. DGB {fr (0O)}f transforms into a DOGB 
{ff (t O J, O)}^, the energy of which rapidly grows with 
frequency due to the tails. 

Increasing the coupling, the dependence H{Auj) can 
become non-monotonous (e.g. C = 0.17 in Fig. [SJ, which 
is connected to an additional transition from discrete- to 
continuous-like breather solution p^l- Close to the AC 
limit (C ^ 1 or 3> wi), the lower-band sub-field of 



the breather is well localized, and the two peaks of the 
upper-band sub-field (to the left and to the right of the 
breather center, see Fig.|21l are situated on the 'light' sites 
nearest to the breather center. We call this configuration 
'discrete-like' breather. By contrast, close to the contin- 
uous limit the lower-band sub-field extends over many 
lattice sites, and the peaks of the upper band sub-field 
are located far from the breather center - 'continuous-like' 
breather. Consequently, the continuation of a breather 
from the anti-continuous to the continuous limit is always 
accompanied with an infinite number of tranformations, 
where each transformation corresponds to a 'jump' of up- 
per band sub-field peaks from positions n = Nc^tiq to po- 
sitions n — Nc±{no + 2) with no = 1, 3, 5, ... for symmetric 
breathers and uq = 2, 4, 6, ... for antisymmetric breathers. 

In Fig.Eldifferent regimes of transitions from 'discrete- 
like' breathers with peaks of the upper band sub-field at 
n = Nc ± 1 (see insets in Fig. |2Ia),(d)) to 'continuous- 
like' breathers with peaks of the upper band sub-field at 
n — iVc±3 (insets in Fig.|nib),(e)) are demonstrated in the 
plane {H, Auj) for continuation in frequency at different 
fixed values of coupling. For small enough C, a 'discrete- 
like' breather smoothly transforms into a 'continuous-like' 
one while decreasing the frequency towards the lower gap 
boundary {Auj —2(5^). At a certain critical coupling 
Ci the derivative dH/duji, turns to zero at the transition 
point (see Fig. EJa)). Further increase of the coupling will 
lead to the appearence of a part of the H{Auj) curve with 
a negative slope in the transition region (see Fig.Efb)). As 
a consequence, discrete-like and continuous-like breathers 
will possess real instabilities of the Vakhitov-Kolokolov 
type l^T\ near the transition point. These instabilities are 
produced by the collisions of the N2 localized eigenmode 
with its complex conjugate (see Sec. I5.1|l . 

At the second critical value C2 > Ci the slope of the 
'unstable' part of the curve H{Auj) will become vertical, 
and for C > C2 it bends over and the dependence H{Auj) 
becomes multi-valued in the transition region (Fig. |Bfc) 
and lower curve in Fig. Eld)). In this regime the transi- 
tion between discrete-like and continuous-like breathers is 
possible only via an intermediate solution. 

For large enough coupling, multi-breather solutions get 
involved in the transition from discrete-like to continuous- 
like breathers. In the upper right part of Fig. IHId) the 
bifurcation loop for the breathers {it (i O T 0)}s 
(lower part of the loop) and {trjJI (i O t 0)}s (upper 
part) is shown. Being constructed at frequencies above 
the linear spectrum, these DOGBs were continued in fre- 
quency down to the values inside the gap^, where they 
bifurcate with each other close to the bifurcation point of 
the continuous-like DGB with the intermediate breather. 
When the coupling reaches a third critical value C3, these 
two bifurcation points coincide (Fig.|B|[e)), and the four bi- 
furcating solutions are identical in that point. For C > C3 
the bifurcation picture drastically changes. The curves 
for a continuous-hke DGB {ff (0O)}f and for a DOGB 

® Note that these two solutions cannot be continued in fre- 
quency inside the gap at small values of coupling. Therefore, it 
is impossible to associate any 'gap' coding sequences to them. 
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Fig. 6. Transitions from 'discrete-like' (larger AiS) to 'continuous-like' (smaller Auj) breathers. The dependencies of breather 
Hamiltonian energies Q on frequency detunings Auj at different fixed values of the coupling: (a) C = 0.169; (b) C = 0.18; 
(c) C = 0.185; (d) C = 0.191; (e) C = 0.192; (f) C = 0.28. The insets in (a),(b),(d) and (e) show the structure of 'discrete-like' 
DGB (C = 0.169, Zio; = -0.05), 'continuous-like' DGB (C = 0.18, Zio; = -0.1), 'discrete-like' DOGB (C = 0.191, Zioj = 0.01, 
lower curve) and 'continuous-like' DOGB (C — 0.192, Ald = 0.01, middle curve), respectively. In all cases, 5^ — 0.1 and A'^ — 242. 



O 



{1T J| (I O t 0)15 coincide, so that a single-breather so- 
lution DGB {-fl- (OO)}^ can be smoothly continued up to 
frequencies inside the linear spectrum (upper solid line in 
Fig-EJe)). At the same time the curves for the intermedi- 
ate breather, connecting the continuous- like and discrete- 
like DGBs, and for the DOGB {frilL (i O T O)}^ also 
coincide (dashed line in Fig.Efe)). Therefore, the discrctc- 
hke DGB {fr (0O)}f connected to the DOGB {ff (T O i 
0)}g now bifurcates with the DOGB (j O t 0)}g 

at frequencies inside the gap (lower solid line in Fig.|^e)). 

Further increase of C will lead to the movement of the 
bifurcation point of the discrete- like DGB {ff (OO)}^ and 
the DOGB (i O T 0)}% towards the upper bound- 

ary of the gap. Finally, above the forth critical value of the 
coupling C > C4, the discrete-like DOGB {fT (T O i 
as weU as the DOGB {^i^ (i O T 0)}g cannot be contin 
ued in frequency inside the gap, since they bifurcate with 
each other above the upper boundary of the gap Ald = 
(see Fig.EJf)). Consequently, for such large couphng only 
continuous-like DGBs {ff (OO)}^ can exist in the gap. 

Similar transitions from discrete- to continuous-like 
breather solutions also occur for the continuation in cou- 
pling of DGBs and DOGBs at a fixed value of the fre- 
quency detuning Alo (I19|I . In the lower half of the gap 
(—2(5^ < Auj < —6'^) the transitions are smooth. In the 
upper half of the gap the energy curve H{C) bends over in 
a similar way as in Fig.l^Ic), and the discrete- like breather 
transforms into a continuous-like one via an intermediate 
solution. Above the upper gap boundary {Alo > 0) the 
DOGBs ^ (i O T 0)}g and {^[^ (i O T 0)}° get in- 
volved in the transition between discrete- and continuous- 



like solutions (Fig. . In the lower left part of Fig. Ufa) 
the bifurcation loop for these two solutions is shown, while 
the upper right part shows the loop for transition from 
the continuous-like (upper solid line) into the discrete-like 
(lower solid line) DOGBs {ft (T O i 0)}g via the in- 
termediate out-gap breather (dashed line). At a certain 
critical Alo the bifurcation point of the DOGBs {fT JJ- (i 
O t 0)}s and {triJI (i O T coincide with that of 

the continuous- like and intermediate out-gap breathers. 
At higher frequencies the continuous-like DOGB can be 
continued smoothly in coupling down to C = 0, with the 
coding sequence {fl" -IJ- (| O t O')}^ (upper solid line in 
Fig.I^b)), while the discrete-like DOGB {tf (T O i 0)}^ 
bifurcates with the DOGB {friJJ. (j O t 0)}g (lower solid 
and dashed lines in Fig.l^b), respectively). An analogous 
scenario will then be repeated at higher frequencies, for 
the transition from a solution with peaks of the upper 
band sub-field at n = iVc ± 3 to one with upper band 
sub-field peaks at n = A^c ± 5. In Fig. C{b) one can see a 
bending of the H{C) curve for the continuous-like solution 
at C ~ 0.68, associated with this subsequent transition. 

A similar switching of the bifurcation scenario with 
change of coupling was also observed for 'phantom' brea- 
thers associated with second-harmonic resonances in mono- 
atomic KG chains |39| . 



5.3 Bifurcations of DOGBs with different tails 

As described above, gap breathers {ff (OO)}^ gain non- 
decaying tails with wave number q = ix jl when continued 
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Fig. 7. Transitions from 'discrete-like' to 'continuous- like' out-gap breathers (see text for explanations). The dependencies of 
breather Hamiltonian energies H ^ on the coupling C at different fixed values of the frequency detuning Aco: (a) Aco = 0.1 • (5^ 
and (b) Acu = 5^ . The values of the system parameters are the same as in Fig. |S| 




Fig. 8. Hamiltonian H lO vs. frequency uib for different solutions. Vertical lines show the boundaries of the upper branch of 
the linear spectrum lu2 = 1.14, tuu ~ 1.148. Curves in (b) show, from top to bottom: continuation DGB {-fi (0O)}g DOGB 
{lt (T O i O)}^, bifurcation loops for DOGBs with 1,2,4,8 and 12 'defects' (see text for explanation). 5^ = 0.1, C = 0.02, N = 50. 



into 'normal' DOGBs above the upper gap boundary, and 
their energies rapidly grow. However, there exist also local- 
ized breathers with frequencies above the linear spectrum 
- 'on-top' breathers (DOTB), with exponentially decaying 
tails of wave number q = n (for stationary solutions) . The 
simplest example is the single-site DOTB with coding se- 
quence {tr (0O)}g. Attempting a continuation of this solu- 
tion in frequency down to the linear spectrum, it bifurcates 
with another on-top breather {^l (OO)}^ rather far from 
the 'top' frequency uJu of Hnear waves (Fig. (Hfa), lowest 
curves). Another DOTB {fft (OO)}^ (three-site breather) 
can be continued much further down in frequency, bifur- 
cating with the DOTB {frt O T (O0)}g close to the point 
u = oju (Fig.lHIa), middle curves). The slope of the curve 
H{u)t) for the DOTB {frt (O0)}g is approximately the 
same as for the DGB {))■ (OO)}^. In that sense this on- 
top breather, rather than the DOTB {ff (OO)}^, can be 



considered as the 'on-top counterpart' of the gap breather 
(however, this is true only for small coupling, see below). 

Besides the 'normal' DOGBs, another type of out-gap 
breather with tail wave number q = 7r/2 exists above the 
gap - the DOGB {(O t O i)}g (and its antisymmetric 
counterpart {0(0 t O 1)}a)- The dependence H{oJi,) for 
the continuation in frequency of this DOGB at fixed cou- 
pling C = 0.02 is shown in Fig. |Sfa) (dash-dotted hue). 
Originating from the linear wave with g = 7r/2 at the up- 
per gap boundary, it represents a kink-like excitation of 
'light' sites ('phase domain wall') accompanied by local- 
ized excitation of 'heavy' sites at non-zero coupling, and 
corresponds to the kink-like out-gap soliton in the contin- 
uous limit (e.g. [7ll29| ). However, above a critical coupling 
Ccr ~ 0.188 the DOGB {(O t O i)}g does not exist at 
Au) — > 0, since being continued in frequency down towards 
the upper gap boundary at fixed coupling C > Ccr, it bi- 
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furcates with the solution (O i O T)}^ at Auj > 0. 
Instead, at C > the solution {JJ. 0(0 [ O t)}s can 
be continued down to the upper gap boundary, corre- 
sponding to the kink-like soliton in this regime. This is 
similar to the switching of the coding sequence for the 
continuous-like 'normal' DOGB from {H (T O i 0)}s to 
{tr -IL (I O T 0)}g discussed in SecO 

Apart from DOGBs with tail wave numbers q — tx jl 
and DOTBs with q = tt, one can also construct DOGBs 
with any other tail structure, if the chosen wave number 
is allowed for the particular system size. As the coding se- 
quences become complicated, we introduce a simplified no- 
tation. Restricting to symmetric breather configurations 
we omit the subscripts 5, and likewise the superscripts 
O will be omitted as all DOGBs are constructed at fre- 
quencies above the gap. For DOGBs with frequencies in- 
side the linear spectrum the upper band sub- field ('light' 
amplitudes) should have non-zero tails because of reso- 
nances with linear waves. Therefore, the codes for 'light' 
amplitudes should be, generally, non-zero. Restricting to 
DOGBs which are 'single-site' in the lower band sub-field, 
all codes for 'heavy' sites are zero except for the central 
excited site, and the codes for 'heavy' amplitudes can be 
omitted. For example, the simplified coding sequence for a 
DOGB {tr (0O)}g is [(0)], and for DOGB {ff (T O i 0)}g 
is [(ti)] (we use different types of brackets to distinguish 
normal and simplified coding sequences). 

To understand the total bifurcation picture of DOGBs 
with different tails, we construct new solutions by insert- 
ing 'defects' into the tail structure of a DOGB [(ti)]- Con- 
sidering this sequence as the 'right' one, we replace the 
'right' code at some particular place k with a 'wrong' one 
together with introducing an additional phase shift, i.e. 
a change of the signs of all the subsequent codes (i.e. we 
insert a defect into the 'right' coding sequence to obtain a 
standing wave with wave number smaller than q — 7r/2). 
For instance, denoting the number of codes in the sim- 
plified coding sequence as Nl (i.e. the number of 'light' 
sites to the right of the breather center, e.g. — 13 
if = 50), the solution with one 'defect' at the posi- 
tion k (in the simplified coding sequence, corresponding 
to site n — Nc + 'i.k — 1) has the simplified coding se- 
quence: [[ri)i=i^2,...,k-i,Wk,{,-ri)i=k+i,k+2,...,NL\- Here 
is the 'right' code at the i-th position, and Wk is a 'wrong' 
code {wk is '0' or 'J,' if rk is 'f, Wk is '0' or 'f if is 
'J,'). A DOGB [0(ti)] is an example of a solution with a 
defect at the first position, while a DOGB [TiTT (it)] has 
a defect at the 4th postition. Similarly DOGBs with more 
defects can be constructed. A single-site DOTB [(0)] can 
be considered as a solution with Nl 'defects'. 

We find, that being continued down in frequency to- 
wards the linear spectrum, each solution with 'defects' bi- 
furcates with a solution having the same number of 'de- 
fects', with coding sequences differing only at the 'defects' 
nearest to the center. For example, the DOGB [ttit ii 
(Ti)] with three 'defects' at 2nd, 5th and 7th positions 
bifurcates with the DOGB [T it ii (Ti)], with a dif- 
ferent type of 'defect' at the 2nd position. The resulting 
bifurcation picture is shown in Fig. IHl^b) , where bifurca- 
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Fig. 9. The intrinsic structure of the bifurcation loop for the 
DOGB with one 'defect' in Fig.|HIb). Dashed black, solid dark- 
grey, solid light-grey and solid black lines represent solutions 
with one 'defect' at 1st, 2nd, 3rd and 9th positions respectively. 



tion loops for solutions with 1,2,4,8 and 12 'defects' are 
plotted in different greyscales, while the solid black curve 
corresponds to the transition DGB {ff (0O)}f DOGB 
[(ti)]. Note that in Fig. |H1A^l = 13, and the bifurcation 
loop for solutions with 12 'defects' in Fig.lSfb) is the same 
as that for the DOTBs [t (0)] and [tt (0)] in Fig.lH^a). 

Each bifurcation loop corresponding to a particular 
number of 'defects' has a complex intrinsic structure, as 
these 'defects' can be at different positions. For each num- 
ber of defects one can find some 'optimal' configuration, so 
that the solution with 'defects' placed on these 'optimal' 
positions can be continued further down in frequency than 
all other solutions with the same number of 'defects'. The 
loops shown in Fig. |Hl[b) represent such optimal configu- 
rations. In Fig. I^lthe intrinsic structure of the bifurcation 
loop for solutions with one 'defect' at 1st, 2nd, 3rd and 
9th ('optimal') positions, respectively, is shown. 

Increasing the coupling, the bifurcation picture remains 
qualitatively the same. However, the slope of the curve 
H{ujb) for the DGB {fr (0O)}f chan ges, so that for larger 
coupling, when the DGB has 'continuous-like' structure 
with upper-band sub- field peaks at n = Nc ±3, its 'on- 
top' counterpart (with approximately the same slope for 
7J(wb)) is the DOTB {frt O i (O0)}g. Similar changes 
of the DGB's 'on-top' counterpart are expected for a fur- 
ther increase of coupling, as the peaks of the upper-band 
sub-field will move further away from the breather center. 

The bifurcation scheme for out-gap breathers in Fig. 
|Hl[b) is qualitatively similar to that for phonobreathers 
with different tails resulting from second-harmonic res- 
onances in monoatomic KG chains "W. However, in the 
modulated DNLS model out-gap breather tails apparently 
always have non-zero asymptotes, and thus we found no 
analogue to 'phantom' breathers with vanishing tails. 
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Fig. 10. Imaginary parts of eigenvalues iuje of Ijl2|l for symmetric discrete-like (a), intermediate (b), continuous-like (c) and 
antisymmetric (d) 'normal' DOGBs at frequency deviation 1191 1 Au = 0.1- 5^. Light-grey crosses (dark-grey triangles) correspond 
to eigenvalues with k = +1 {k = —1), and black stars to unstable eigenvalues, (a)-(c) correspond to upper right loop in Fig.|7[a). 



5.4 Stability of out-gap breathers 



values with opposite Krein signatures: 



We analyze the stabihty of DOGBs by continuation in 
couphng at fixed frequency detuning (|19|l , focusing on the 
three pairs of most important solutions: 'normal' DOGBs 
{lt (T O i 0)}g, {0 tr (T O i 0)}% being the smooth 
continuation of symmetric and antisymmetric DGBs in 
frequency above the gap for small coupling; DOTBs {tl" 
(0O)}g, {0 tr (0O)}a, being the fundamental Wtop' 



{Auj + 2S^y + 8C^ cos2((7) 



± 



(23) 



breathers; and DOTBs {frt (O0)}g, {0 frt (OO)}^, being 
the 'on-top counterparts' of DGBs for small coupling. We 
fix the system size N = 242 and gap parameter 5^ = 0.1. 



5.4.1 'Normal' DOGBs 

In Fig. llUl we show stability results for symmetric discrete- 
like, intermediate, continuous-like, and antisymmetric out- 
gap breathers at frequency detuning Aui = 0.1 • 6^. (Note, 
that in 'normal' DOGBs the tail wave-number is q = 7r/2, 
so according to ((T^ Auj ~ ujh~uJo{Tr/2) = lj;, — cj2-) As for 
DGBs (Sec. I5.l|l . there are two bands of extended eigen- 



±^^iAuj + 2^2)4 + 16C2 cos2(q) (454 - z\a;2)| ' ^ 

with all wavenumbers q allowed for the particular system 
size. But now the band {iw+j with positive Krein signa- 
ture always has its lower boundary at zero. Moreover, as 
for non-localized solutions Krein instabilities correspond- 
ing to extended eigenmodes survive also for infinite sys- 
tems (e.g.jS?!), overlapping of the extended bands is more 
essential for extended DOGBs than for localized DGBs. 

Increasing the coupling, several localized modes bifur- 
cate from the top of the band with k — +1 {P1-P5 in 
Fig. ^1 (a),(d)) and from the bottom of the band with 
K = -1 (iVl and N2 in Fig. [TUI (a),(d)). Penetrating the 
bands with opposite Krein signature, these modes produce 
oscillatory instabilities. The A^l mode also causes a real in- 
stability with an eigenvector of opposite symmetry to the 
DOGB, colliding with its complex conjugate at zero (PN 
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Fig. 11. Blow-up of a regime close to the origin of Fig. llUf d). 



instability yielding 'exchange of stability^' between sym- 
metric and antisymmetric DOGBs, cf. Sec. 15. 1(1 . Similarly, 
the N2 mode causes a real instability with an eigenvec- 
tor of the same symmetry as the DOGB, in the regime of 
transition from discrete-like to continuous-like solutions. 

However, the behaviour of the localized modes at small 
coupling is different for gap and out-gap breathers. For 
symmetric DGBs and DOGBs, in the limit of small C 
only two modes can stay localized: PI and P2, associated 
with symmetric and antisymmetric oscillations on the sites 
closest to the breather center. In the case of DOGBs, in 
contrast to DGBs, there is no threshold Aujcr, and these 
two eigenvalues bifurcate from the extended band immedi- 
ately when C becomes non-zero, for all DOGB frequencies 
(see Appendix 1X1 Tabled]). The differences between each 
of the PI and P2 eigenvalues and the top boundary of 
the extended band are in this case proportional to \/C at 
small values of C (see Fig. ^| and Table ^ compare with 
C^-dependence for DGBs as discussed in Sec. I5.1|l 

For antisymmetric 'normal' DOGB {0 H (t O i O)}^ 
an additional mode iVO appears (see Fig.lll|l. localized on 
the 'central' site Nc which is the only 'light' site with code 
'0'. Note, that in an antisymmetric DGB {0 (OO)}^ aU 
'light' sites have zero codes, therefore there is no such 
localized mode for breather frequencies inside the gap. 

Another important difference between antisymmetric 
'normal' DOGBs and DGBs is, that while the latter al- 
ways have a real instability at small coupling (see Sec. 
I5.1|l . the former are always stable in a small interval of C: 
[0, Co] (see Fig. 111(1 . A localized mode PO, with eigenvector 
of opposite symmetry to the DOGB, now bifurcates from 
the upper boundary of the positive Krein signature band, 
similarly to the PI and P2 modes. The eigenfrequency is 
proportional to C for small C (see Table^in Appendix IXI 
eigenmode {P0,A,O)), and the eigenvector for C ^ in- 
volves oscillations not only in the 'excited' 'heavy' sites, 
but also in the subsequent 'light' sites. At a certain non- 
zero value of C it collides with the iVO mode producing 



an oscillatory instability, associated with a resonance be- 
tween oscillations in the central 'light' site and symmetric 
coupled oscillations in the two neighboring 'heavy' sites. 

The A^l localized mode, which is responsible for the 
'exchange of stability' between antisymmetric and sym- 
metric 'normal' DOGBs, now appears only at large cou- 
pling, in the continuous-like regime. 



5.4.2 'On-top' DOGBs 

As indicated in Sec. 15.31 'on-top' breathers cannot be con- 
tinued in frequency down to the linear spectrum^. Anal- 
ogously they cannot be continued in coupling up to ar- 
bitrary values at fixed frequency detuning Au, since they 
bifurcate with other DOTBs in the same way as when con- 
tinued in frequency at fixed C (see Fig.lHJa)). By contrast, 
the DOTBs can be smoothly continued by simultaneously 
increasing the frequency detuning and the coupling at 
fixed value of the ratio Auj/C. This corresponds to the 
continuation towards the 'non-modulated' limit S'^/C 
(see Eq. dJ), i.e. towards usual 'on-top' breathers in ho- 
mogeneous media. Only for = the smooth contin- 
uation in coupling at fixed frequency detuning can be per- 
formed; at any non-zero 6'^/C the exact limit Auj/C 
cannot be reached because of the bifurcations. 

In Fig.^lresults of continuation in coupling of DOTBs 
{il (0O)}O {0 ii (0O)}O, m (O0)}O and {0 M (OO)}^ 
at fixed frequency detuning are shown. Each DOTB is 
continued from C = up to the bifurcation point, where 
it bifurcates with DOTBs {ffi (O0)}g, {0 fri (OO)}^, 



{trt O t (O0)}g and {0 trt O t (OO)}^, respectively. 
For all DOTBs the extended eigenvalues have the same 
(negative) Krein signature. Still they can be divided into 
two bands ('upper' and 'lower'): 



iuJlo =UJb- t^a(g), 



(24) 



Krein instabilities generally also exist here, see Fig. 1101 



where Wo,a(9) are frequencies from the upper and lower 
bands of the linear spectrum 

Close to the bifurcation point, all the described DOTBs 
possess two real instabilities with symmetric and antisym- 
metric eigenvectors, respectively. These instabilities are 
similar to those produced by the A^l and A^2 modes in 
the case of DGBs and 'normal' DOGBs in the region of 
transition from discrete- to continuous-like solutions. 

Like 'normal' antisymmetric DOGBs, antisymmetric 
DOTBs have no instabilities at small coupling, so for each 
pair of symmetric and antisymmetric DOTBs there is an 
interval [0, Co] of simultaneous stability. The localized mode 
associated with symmetric coupled oscillations in the two 
'excited' 'heavy' sites (similar to the PO mode for a 'nor- 
mal' antisymmetric DOGB) bifurcates from zero and, hav- 
ing K > 0, collides with another localized mode associ- 
ated with oscillations in the 'central' 'light' site (similar 

* We here only discuss DOTBs with nonzero oscillations at 
'heavy' sites as C — » 0. There are also solutions with nonzero 
codes only at 'light' sites, e.g. {t (OO)}^, which can be con- 
tinued down to the linear spectrum. 
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Fig. 12. Imaginary parts of eigenvalues iue. of l|12|l for continuation in coupling of DOTBs {-(^ (0O)}g (a), {0 "1^ (0O)}2 (b), 
{1>T (O0)}g (c), {0 I^T (O0)}2 (d) at frequency deviation (O Auj = 5^ (note, that g = tt for DOTBs). Light-grey crosses 
(dark-grey triangles) correspond to eigenvalues with ft = (k = —1), and black stars correspond to unstable eigenvalues. 



to the iVO mode) . This resuhs in an osciUatory instability 
of antisymmetric DOTBs (Fig. IT^b). (d)). The DOTBs 
{trt (O0)}g and {0 frt (O0)}2 have two additional local- 
ized modes with k > 0, similar to the PI and P2 modes 
described in Sec. 15.4.11 Penetrating the extended bands 
or colliding with other localized modes (with k < 0) they 
also produce oscillatory instabilities (Fig. I12r c). (d)). 



6 Conclusions 



not possess any oscillatory instability when the breather 
frequency is in the lower half of the gap. 

The transition from gap to out-gap breathers was stud- 
ied by varying the breather frequency at different fixed 
values of the coupling. Such transitions were found to oc- 
cur in different ways, depending on the value of the cou- 
pling, and a relation to the transition from discrete- to 
continuous-like breathers (described earlier for DGBs in 
the diatomic KG model |Ej) was revealed. 



We have investigated the linear stability properties of dis- 
crete gap and out-gap breathers in the modulated DNLS 
model, which in particular describes electromagnetic wave 
propagation in an array of weakly coupled optical waveg- 
uides of different widths. The coupling constant was varied 
from C = (anti-continuous limit) up to values, at which 
the localization length of the breather becomes much larger 
than the array constant (continuous limit). 

Different types of oscillatory and real instabilities of 
DGBs were revealed, similar to those described for DGBs 
in the diatomic Klein-Gordon model ^7]. However, in con- 
trast to the KG case, DGBs in infinite DNLS chains do 



Discrete out-gap breathers with tails of different wave 
numbers were investigated, and a general bifurcation pic- 
ture described. DOGBs were shown to, generally, pos- 
sess the same types of oscillatory and real instabilities 
as DGBs. Still, there are important differences between 
stability properties of DGBs and DOGBs. Unlike anti- 
symmetric DGBs, antisymmetric DOGBs are linearly sta- 
ble for small coupling, so there is an interval in coupling 
[0, Co] where symmetric and antisymmetric DOGBs are 
stable simultaneously. Also, since DOGBs generally are 
non-localized, they possess strong oscillatory instabilities 
associated with resonances between extended lattice modes 
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of the upper- and lower-band sub-fields. For DGBs such 
instabilities do not survive in the limit of infinite chain. 

The process of breather frequency penetration into the 
linear spectrum and bifurcation of a DGB into a DOGB 
was found to have many qualitative similarities with the 
breather second harmonic penetration into the linear spec- 
trum and formation of 'phonobreathers' with non-zero 
tails, discussed for a monoatomic KG chain in [3^1 • How- 
ever, we found no analogue to 'phantom' breathers with 
vanishing tails, since all DOGBs with tails of different 
wave numbers apparently arc non-localized at frequencies 
inside the linear spectrum. The question about the exis- 
tence of such localized excitations inside the linear wave 
spectrum in modulated models is thus still unresolved. 

Another interesting and important question, which we 
address for future investigations, concerns the possibility 
of finding exact moving gap and out-gap breathers. As we 
have shown here, DGBs and DOGBs in the modulated 
DNLS model exhibit the 'exchange of stability' between 
symmetric and antisymmetric modes, and thus a consid- 
erably increased mobility of these excitations is expected 
for parameter values close to the stability exchange point 
We hope that our results will stimulate experimen- 
tal activity to directly observe this phenomenon, e.g. in 
coupled waveguide arrays. 

We recently demonstrated such DGB mobility in the 
diatomic KG chain 17 , and connected it with a vanishing 
of the PN barrier. However, the determination of the PN 
potential, and the PN barrier as the energy difference be- 
tween symmetric and antisymmetric breathers, is not un- 
ambiguous in the KG model, since for fixed values of the 
system parameters there exist continuous families of sym- 
metric and antisymmetric DBs with different frequencies. 
In the modulated DNLS model, the existence of an addi- 
tional conserved quantity, the excitation norm, makes it 
natural to define the PN barrier as the difference between 
energies of symmetric and antisymmetric DBs with the 
same norm^Jj. Thus it is of interest to study more in de- 
tails the 'exchange of stability' in the modulated DNLS 
model, and to look for exact moving breathers in this 
regime. 

Another interesting effect, which also should be ex- 
perimentally observable, arises as a consequence of the 
existence of DGB internal modes with frequencies above 
the continuous spectrum (e.g. P1-P5 in Fig.|3|). Since then 
also all multiples of the mode frequency will lie outside the 
continuous spectrum, similar arguments as in |43) prove 
the existence of exact localized quasiperiodic gap breathers, 
with two, generally incommensurate, frequencies. These 
solutions are characterized by persistent large-amplitude 
internal intensity oscillations. It is particularly interest- 
ing to note, that since e.g. the PI mode of the symmetric 
DGB is always above the continuous spectrum in the lower 
half of the gap when —1.56 < Auj/ 6'^ < —1, quasiperi- 
odic solutions bifurcating from this mode should exist in 
continuous families, arbitrarily close as well to the anti- 
continuous as to the continuous limit. The properties of 
such solutions, which cannot exist in KG models, will be 
discussed elsewhere. 
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A Calculation of localized modes at small C 

For localized eigenmodes at small C, the eigenvalue prob- 
lem H15|l can be simplified. Indeed, at small coupling these 
eigenmodes are mainly localized on a few central sites 
(if the corresponding eigenvalues lie inside the extended 
band, these eigenmodes are not purely localized, but the 
central site amplitudes are much larger than the tails). In 
this case one can put {a,b)n = ior \n — Nc\ > N^, where 
iVf characterizes the localization length of an eigenmode. 
Moreover, using C as a small parameter, one can obtain 
approximate expressions for the central site amplitudes of 
the unperturbed solution for DGB {ff (0O)}f : 

02^^ = 26^ + Auj- -^C^ + O (C^) , (25) 
DGB {Otr (0O)}${: 

4'%^ = 0, 

±1 = 25^ + Alj~ -^C^ + O (C^) , (26) 
Auj 

o 2(5^ + Auj o , 

DOGB U (T O i 0)}g: 

0k = 26' + AUJ + 2^11^^ + O (C^) , (27) 

and DOGB {0 (t O i O)}^: 

0k±i = 26' + ^^+ \l 26^^+Auj ^ + ^ ' 

Note that Auj < {Auj > 0) for gap (out-gap) breathers. 

Thus, in the case of a symmetric DGB {fl" (0O)}g, 
the PI mode is mainly localized on the central 'heavy' 
site with n — Nc, the two neighboring 'light' sites with 
n = Nc ± 1, and the subsequent pair of 'heavy' sites with 
n — Nci2. Therefore, for small C one can put (a, 6)„ = 
when jn — A^d > 3 in equations (|15|l . Substituting the 
approximate breather solution l|25(l into (|15|l and using 
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Table 1. Localized eigenmodes at small coupling. Eigenfrequencies are given at lowest significant order of C, with the notations: 
Ai = 23^ + Ah) (i.e. frequency detuning measured from the lower gap edge), = |^165^ + Zif ± \J 1285* + Zif j / [4Z\(jj^(5^] . 
Schematic eigenvector structures are shown for the central part; higher arrows denote sites with nonzero oscillations at C ^ 0. 



Eigenmode 


Schematic eigenvector structure 


Eigcnfrequency (j;^ 


AlOct 


(PI, 5-, G) 


...oo-ff T ^ T frOO... 


-zio; + 2G^ (4(5^ + Zio;) j Au^ 


-1.56 -(5^ 


(P2,5,G) 


...OOfr i O T liOO... 


-Alo + G" (4(5'' + Ai) 1 2S'Alu'' 


-1.1 ■ 


(P1,A,G) 


...004 T 4 T 1). T ^00... 


-Au + G^^i(+) 


-1.85 • 5'^ 


(P2, A,G) 


...OO4 T i ^00... 


-Zio; + G" (12(5^ + 4zia;(5^ + zio;") / 2Auj^S^ 


-1.44 • 5'^ 


(P3, A,G) 


...oo^t T -It T 4 T ^00... 


-Auj + G^^i(_) 


-0.86 ■ (5'' 


{N\,A,G) 


...OOt 1t T 1> tOO... 


i2Cy/~Al/-^-Auj 


-2-5' 


(Pl,5,0) 


...004 iin j>oo... 


{Q^JAujA^cY" 





(P2,5,0) 


...OOfr i T J>0O... 


{2y/AuoAiCY"-'' 





(P1,A,0) 


...00 T4 0i>i 00... 


{'i^JAujA^cy-" 





(P2,A,0) 


...OOffil>Tl>ifrOO... 


{'i^AujA^cY" 





(P0,A,O) 


...OO4 Ti> T 1>T 4OO... 


C^J2A\ - Aw■^|^JA^^Ax 







...OOit t -friOO... 


Alo + 0{C*) 






the symmetry of the PI mode: {a,b)N^-i = (a, 6)Ar^+i, we 
obtain the corresponding eigenfrequency as 

^{Pl,S,G} ^ + ^ (4^2 + + (^3) _ (29) 

Here the superscript marks the type of mode (PI or P2), 
the symmetry of the breather </)i°^ {S [A] for symmetric 
[antisymmetric] breathers), and the type of breather (G 
and O for gap and out-gap breathers, respectively). 

Analogously, for the antisymmetric P2 mode we have: 



{P2,S.,G} ^ 



-Alu 



(30) 



To see if a particular eigenmode is localized or not at 
small C, i.e. whether the corresponding eigenvalue lies in- 
side or outside the extended band, we check the difference 
between the eigenvalue and the band edge ^^j^'^^^^'i 

^{ma.,G} ^ ^^(^^ _ ^ _ ^2 ^ ^^4 + 4^2^ (31) 

from which it bifurcates (see Fig. |3J|. This yields for the 
P\ mode 



,{P1,S,G} 



{maa;,G} 



(32) 



Auj'^6'^ 



Consequently, when the frequency detuning of a symmet- 
ric DGB is above the critical value Auj > Aulr ' ~ 
— 1.56 • (5^, the PI eigenmode is localized at any non-zero 

value of C. By contrast, when Auj < Aucr ' , the PI 
eigenmode is not localized at small values of C. 
Similar results are obtained for the P2 mode: 



{P2,S,G} 



,G} 



(33) 



C2 



2Auj^S^ 



iS^Auj - 3Alj^) + O (C^ 



The critical value of the frequency detuning for this mode 
is equal to Aujcr^'^'^^ w — 1.1 • iJ'^. 

Using the same technique, expressions for the eigen- 
values and eigenvectors corresponding to localized modes 
at small coupling can be obtained for antisymmetric DGB 
{0 tr (OO)}^, and for out-gap breathers {fT (t O i 0)}g, 
{0 li (T O i O)}^. The results are presented in Table □ 
where for each eigenmode the schematic eigenvector and 
eigenvalue are given at the lowest significant order of C, 
together with the approximate value of the critical fre- 
quency detuning Aujcr-, above which the mode is local- 
ized at small C. Note in particular that the antisymmet- 
ric DGB {0 11 (OO)}^ has an imaginary eigenfrequency 
corresponding to the A^l mode, and thus it is unstable for 
small C. All other eigenfrequencies are real, proving the 
linear stability of the other three solutions for small C. 
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